Generation of arbitrary symmetric entangled states with conditional linear optical 
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An approach for generating the entangled photonic states |^i, ^2) ± |^2, ^1} from two arbitrary 
states and ^2} is proposed. The protocol is implemented by the conditionally induced beam- 
splitter coupling which leads to the selective swapping between two photonic modes. Such coupling 
occurs in a quantum system prepared in the superposition of two ground states with only one of 
them being involved in the swapping. All the entangled states in the category, such as entangled 
pairs of coherent states or Fock states (NOON states), can be efficiently produced in the same way 
by this method. 

PACS numbers: 42.50. Dv, 42.50.Pq, 42.50.Gy, 03.67.Lx 



Bipartite symmetric entangled states refer to a generic 
type in the form \^>i, ^/ 2 ) ± \^2, *i) up to a normaliza- 
tion factor. Such entangled states include the symmetric 
entangled coherent states (SECSs) la, /3) ± \/3, a) [l[ and 
the NOON states \N, 0) ± |0, N) HQ. Both of them have 
found important applications in quantum metrology; see, 
e.g. [J, [5[. A SECS of light fields can be transformed to 
a photonic Schrodinger cat state I7) ± | — 7) .6] simply 
by a beam-splitter (BS) operation. Cat states of mat- 
ter wave and even light field have been experimentally 
demonstrated [7-9], but a photonic one with the suffi- 
ciently large size I7I is still beyond the reach. 

Since the seminal work of Yurke and Stoler [l(| , the ap- 
plication of Kerr nonlinearity has been suggested as the 
direct way to entangle light fields or construct photonic 
cat states [l| . Realizing strong coupling between photons 
via the suitable nonlinear media is, however, a rather dif- 
ficult task. This barrier stimulates the parallel researches 
on creating the approximate states by squeezing (see, e.g. 
[nLllIi and the reference of @) and exploring the proper 
use of weak Kerr nonlinearity (see, e.g. [la - [l5j V 

A less noticed problem with Kerr nonlinearity and 
squeezing is the availability of their single-mode ver- 
sions, which are the basis for all relevant schemes thus 
far. A realistic photonic pulse carries multiple modes 
represented by the field operator £(z) = J\^fe e (in 
one-dimensional space for illustration). For instance, 
under the action of a multi-mode self-Kerr Hamilto- 
nian J dz[£ ^{z)£ (z)) of the unit coupling constant or 

its equivalent form J2ki,k2.k3 ^li-fc2+fe3^i"L"fc3 m tnc 
wave-vector space, the output states can be significantly 
different from the proper ones that should have evolved 
under the sum of single- mode actions ^ k ifl]fik) 2 , even 
if the inputs are exactly single-mode ones. This effect 
of mode entanglement or mode mixing has been studied 
in [lfl [l7[. A consequence of the effect is a very lim- 
ited clean cross phase (similar to that obtained from the 
single-mode cross-Kerr model) under highly demanding 



conditions [181 ] . On the other hand, a multi-mode squeez- 
ing action YlkiHfl-k + ®kQ>-k) of one field as well devi- 
ates from its single-mode version. In contrast, the multi- 
mode BS Hamiltonian Hbs for two fields £\(z) and £2(2) 
takes the form J dz{£l£ 2 (z) + £i£\{z)} = ^2 k {a\ k a,2,k + 

a\ t k<i\ k)i a sum of the individual mode actions. This 
BS coupling enables a multi-mode photonic state to be 
transformed ideally like a single-mode one, because the 
decomposable evolution operator exp{— UHbs} with re- 
spect to the wave-vector modes k acts independently on 
each mode. 

In this paper we provide a method for generating arbi- 
trary symmetric entangled states out of light fields based 
only on such clean BS coupling . Unlike the common 
linear optical setup, the BS coupling we need acts condi- 
tionally on the part of a superposition of quantum states 
at the same spatial location. Below we will show how to 
produce a symmetric entangled state with a conditional 
BS coupling and will give an example of the realization of 
the given type of interaction in a proper quantum system. 

Protocol. — From now on, we use the term mode in the 
meaning of a single wave- vector or a frequency mode, 
since we will consider BS type coupling only. The two 
arbitrary states ^2) we will entangle are treated 

as the single- frequency ones. 

To entangle the two states, we also need an ancilla 
quantum system with two stable states \m) and \g). This 
system can be an atom, as well as an ion, a quantum 
dot or a superconducting qubit. The ancilla system is 
initially in the state |ra), setting the initial state for the 
total system as |<j>o) = I ^1 7 ^2) \m). Then we perform 
a a x rotation between the \m) and \g) and transfer the 
system to the superposition 



|$1> = 1*1, *2> 



lfl»/V2. 



(1) 



Such o~ x rotation can be realized by applying a resonant 
it/2 pulse to the transition \m) — ¥ \g). 
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The above superposition of |m) and \g) works as a 
logic control on the swapping between two input pho- 
tonic modes: the swapping between the photon modes is 
activated in the \g) subspace and does not happen in the 
|m) subspace. Such conditional swapping can be realized 
by the BS transformation 



ai(t) 
62 (t) 



cos xot 



ie 1 * sin xot 
cos xot 



ai(0) 
a 2 (0) 



(2) 

for the time t s — 7r/(2%o); where aj is the annihilation 
operator for the i-th mode and xo is the effective BS 
coupling constant. The BS transformation can be imple- 
mented via the dispersive parametric three-wave mixing 
(TWM) or four- wave mixing (FWM) pro- 

cess. The conditional swapping results in the state 

|4> 2 ) = (|*i, *a) H - 1 1* 2 , *i) \g)) /V2. (3) 

Then, again we perform a a x rotation between \m) 
and \g) to have them transformed as \m (g)) — ¥ \m (<?)) — 
i \g (to)), leading to the following state 



|* 3 )=2[(|*i,*a>-|*2,*i»N 
-i(|*i,*a) + |*a,*i))|s). 



(4) 



Finally, by measuring \m) and \g) (see the method in 
the following example), we make the photonic sector of 
the total state collapse to the target symmetric entangled 
states |*i, * 2 ) ± 1*2, *i)- 

A candidate for the ancilla system should satisfy two 
requirements. First, the quantum system should have 
two long-lived and well separated states between which a 
a x rotation can be performed. The second requirement is 
specified by the swapping stage — the system should have 
an appropriate energy level structure for the formation 
of the TWM or FWM interaction loop where two of the 
transitions have to be strongly coupled to the input fields. 
These conditions can be satisfied by certain trapped nat- 
ural atoms or ions, single color centers, quantum dots 
or superconducting qubits based on the Josephson junc- 
tions, which have multi-level structures and can also be 
strongly coupled to the suitable field modes. 

The idea to induce the conditional interaction of mat- 
ter wave (ion or atom) state superposition with one opti- 
cal mode was proposed in [23| for creating the cat states 
041) H3 ■ The difference of our setup is that it allows to 
realize a conditional coupling directly between two pho- 
tonic modes for their swapping. This is necessary for con- 
structing a SECS \a, (5) ± |/3, a) with a>^ (for making 
a cat state of large size) or a NOON state. Our method 
aims to generate all such states in a unified way. 

Example. — Here the single ion of calcium 43 Ca + 
trapped in ion trap and embedded in optical resonator is 
considered as the ancilla system. The energy level struc- 
ture of the calcium ion 43 Ca + is illustrated in Fig. 1. 
The research on trapped ions is reviewed in [25[, and a 
recent experiment with the setup similar to our proposal 
is reported in [261 ] . 
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FIG. 1: (color online). 43 Ca+ energy level scheme. 



The two ground states |m) and \g) we use are 
45 1/2 (F = 4, to j = 0) and 3D 5/2 (F = 6, mj = 0), re- 
spectively. These particular levels are chosen as the 
ground states for two reasons. First, both of the states 
are long-lived (up to the order of 1 s). Second, due to 
the selection rule and large energy difference between 
them, one of the ground states is excluded from a para- 
metric FWM loop so that the conditional BS coupling 
can be realized. The ion transition 3_D 5 / 2 (-F 1 = 6) o 
4-P3/2 (F = 5) a t 854 nm is coupled to the a + and a~ 
circularly polarized modes prepared in states |*i) CT + and 
1*2)0— ■ I n case when the input modes are in the coher- 
ent states the SECSs will be generated and for the Fock 
and vacuum inputs the NOON states will be obtained. 
As follows from the interaction configuration presented in 
Fig. 2 only the ground state \g) (3D 5 / 2 (F — 6, raj = 0)) 
becomes coupled to the optical modes but there is 
no coupling to the optical modes for the state \m) 
(4Si/ 2 (F = 4, to j = 0)) . In order to perform a a x ro- 
tation between |to) and \g) and bring the ion into the su- 
perposition state in Eq. ([1]), a resonant 7r/2 laser pulse at 
729 nm with tt - polarization is applied to the quadrupole 
transition AS1/2 O 3D 5 / 2 - During the swapping stage 
there are also two classical pumping pulses with the 
orthogonal circular polarizations applied to the transi- 
tions 3£ 3/2 (F = 5, mj = 0) o 4P 3/2 (F = 5, mj = ±1) 
at 849 nm, while the cavity modes a + and a~ are 
coupled to the transitions 3Z3 5 / 2 (F = 6, mj = 0) «-> 
4-P3/2 {F = 5,mj = ±1) in the parametric FWM loop. 

To realize the parametric BS coupling, all real transi- 
tions in the loop should be suppressed and ideally the ion 
should stay in its ground state during the swapping pro- 
cess. Therefore all fields should be highly detuned from 
the resonance and satisfy certain conditions (see the dis- 
cussion below) . By controlling the duration of the clas- 
sical pulses we can control the precise parametric inter- 
action time for obtaining the state in ([3]). After another 
7r/2 laser pulse at 729 nm with 7r - polarization there will 
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FIG. 2: (color online). Level scheme for realizing the con- 
ditional BS coupling between two photonic modes ai and S2 
(blue arrows). The red arrows represent the classical cou- 
pling fields. The white and black circles on levels \g) and 
\m) indicate that they are in a superposition and only one of 
them is involved into the parametric loop. The energy lev- 
els of this general scheme correspond to those of 43 Ca + for 
our example as follows: \m) — > 4Siy 2 (F = 4, mj = 0); \g) — > 
3As/2 (F = 6, mj = 0); \a) -> 4P 3/2 (F = 5, mj = 1); \b) -> 
4P 3 /2 (F = 5, mj = -1) and |c) ->■ 3.D 3/2 (F = 6, mj = 0). 

be the state in The detection of the ground states for 
the final projection onto the target states is implemented 
by exciting the transition 4Sj m —> 4Pjy 2 at 397 nm; see 
the similar technique in [U S| . The presence of the flu- 
orescence collapses the ion wave function onto |m) and 
the absence of the fluorescence indicates the state \g). 

BS mechanism. — The dispersive FWM process for re- 
alizing the conditional swapping in our protocol can be 
implemented in any system with the level scheme in Fig. 
2. The Hamiltonian for the process shown in Fig. 2 takes 
the form (h = 1) 

H = -Aicr aa - A 2 cr bb - 5<7 CC + gia\(Tga 

+ g2ale~ lAFt a gb + £lia ca + fi 2 cr cb + h.c. (5) 

in a rotation frame. Here cry = |i) (j\ is the atomic 
spin flip operator; gi(2) is the coupling constant; Ai( 2 ) = 
w ai(2) ~ w ga(fc) is the one-photon detuning, 5 = u>i +ui a i — 
Lo cg is the two-photon detuning, and Ap = u)\ + u> a i — 
uj'2 — is the four-photon detuning, with w 1 ( 2 ) being the 
frequency of the classical pumping pulse with the Rabi 
frequency ^1(2), ^ai(a2) being the frequency of the input 
pulse. The Schrodinger equation for each energy level 
component of the state |*(t)} reads: 

ij t (ffl* (*)) = 9ia\ (o|* (t))+g 2 ale'^ t (6|* (t)) (6a) 

*!<W)) = -A a <6|*(t)> 

+n 2 (c|* (t)) + g2a 2 e lAFt (g\9 (*)) (6b) 



i-(a|*(*))=-A 1 (a|*(t)) 

+ ni{c|*(t))+ Sl oi{ S |*(t)) (6c) 

*^<c|*(t))=-*<c|*(t))+n 1 <a|*(t)) 

+ n 2 <6|*(t)). (6d) 

The effective BS Hamiltonian for the similar dispersive 
FWM schemes can be derived by the time-independent 
perturbation method [22| . Here we apply the more gen- 
eral method of the adiabatic elimination [29( to show the 
realization of the effective BS coupling. It is important to 
mention that the one- and two-photon detunings should 
be high enough to prevent any real transition from its 
ground state. It is therefore possible to see the effec- 
tive dynamics of the photonic modes while the system is 
staying in the ground state \g). 

First, assuming that initially the system is prepared 
in its ground state \g), i.e. (&|*(i )) = (a|*(f )) = 
(c|* (to)) = 0, we eliminate the transitions from the state 
\g) to I a) and \b). Under this assumption we integrate 

Eqs. dig) - (unm 

(g\* (*)) = (.gl* (t Q )) - l f dr[ 9l a\ (o|* (r)> 

+ g 2 ale- iA ^{b\*(T))} (7a) 

(6|*(i)) =-* / dre lA2( *- T) 

«/ to 

x [n a (c|* (r)> + 52 S 2e tA ^ ( 5 |vK (r))] (7b) 

(o|*(i)) =-* / dre tAl ^ T » 
•/to 

x [fii(c|*(r))+ 5l ai<fl|*(r))]. (7c) 

Then we substitute Eq. (JTa} into Eqs. (JTHJ) and (JTcJ to 
obtain the relation 

m (t)) = ^ <c|* (f)) + 92 ^ Ft (,9|* (*)> (8a) 

(a\V(t)) = ^{c\*(t)) + ^{g\V(t)) (8b) 

where we use the assumption gi^fnl/Ai <C 1 and keep 
only the first order of the small term, and n, is the aver- 
age photon number of the i-th input mode. 

Next, in order to obtain the decoupled dynamics of the 
effective two- level system of \g) and \c), we substitute 
Eqs. (|8al) and (|8bf into Eqs. (|6a| and (fod)) and obtain 

^-/^ (g\* (t)) = -ie^fi (t) (c|* (t)) (9a) 
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— e 

dt 



(c|#(t)) = -ie^/3(i) (t)>, 



(9b) 



where we have introduced the following functions a c — 
-5 + 0?/A 1 + 0|/A 2) a s = gfa\ai/A 1 +g$ala 2 /A2, and 
/3 = fiiffiOi/Aa + fi 2fl2 a 2 e lAF 7A 2 . 

The dynamics of the states \g) and |c) will be decoupled 
further. Integrating Eq. (|9b[) we get the relation 



<c|*(t)> = i 



Ai 



<$!*(*)>, (10) 



where we keep only the first order of the parameter 
gi^/ni fli/ (SAi) < 1. 

Finally, substituting Eq. (fl"0)) into (|9a|. we obtain the 
decoupled evolution of the state \g) 



i-(g\*(t))=H eff (g\*(t)) 



(11) 



where 

Heff(t) 



gfa\ai 1 / ^igid\ fi 2 fir 2 a| _ lAj! , t 



Ai 



Ax 



A, 



,2^ 



Ai A 2 y A 2 



The conditions gi^/ril/ Ai <C 1, gi-^nfili/ (SAi) <§; 1 lead- 
ing to the above effective dynamics prevent the one- and 
two-photon transitions out of a ground level and can be 
realized by adjusting the system parameters. For our ex- 
ample using 43 Ca + with gxm ~ 10 MHz, it is possible to 
set the Rabi frequencies f2i( 2 ) ~ 1 GHz, the one-photon 
detunings Ai( 2 ) ~ 1 GHz, and the two-photon detun- 
ing 8 ~ 1 GHz, given the average photon numbers up to 
rii = 100. The symbol "~" means the order of the values 
here. The sizes y/rn of the states to be entangled can be 
made larger simply by increasing the detunings. 

Swapping operation. — The unitary evolution operator 
of the time-dependent effective Hamiltonian H e f / in (TT2|) 
can be decomposed as 

Texp|— i J drff e //(r)^ = exp |— iriia\ait — iJ7 2 a 2 a 2 i| 
x T exp | -ixo J e~ lSFT dr a x a\ + h.c.X , (13) 



where T stands for a time-ordered operation, r)i = ^ 



~SAf ' X° — 5AiA 

general form of such decomposition is given in [3 
first of the decomposed operators in Eq.(fT3|) is a phase 
shift operation and the second is a BS operation. For 
example, by tuning the system parameters so that the 
conditions Sp = and S = — 2Ai (assuming g\ = g 2 
and T 3 - = = 1) are satisfied, their combined action 
implements an ideal swapping a\ O a 2 after the time 
t s accumulating Ixo^sl — 0.57T. Given the data following 



= and <5 F = A 



f + m - r i2- 
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FIG. 3: (color online). Upper: (a) Relation between the swap- 
ping time and the ratio xo/Sf- Lower: (b) Fidelity of the 
generated states with the target output ± /3) | a) with 

a = 6\/2 and /3 = y/2, which can be converted to the cat 
states I7) ± I — 7} of 7 = 5. A unit fidelity can be reached in 
the regime |»?i/xo| < 1 at a not so large ratio xo/Sf- 



Eq. (112[) a pair of input states could be entangled within 
a few microseconds. 

(color online) Upper: (a) Relation between the swap- 
ping time and the ratio Xo/^f- Lower: (b) Fidelity of the 
generated states with the target output \ce)\/3) ± \/3)\a) 
with a = 6 y/2 and f3 — \/2, which can be converted to 
the cat states I7) ± | — 7) of 7 = 5. A unit fidelity can be 
reached in the regime |?7i/xo| < 1 at a n ot so large ratio 

Xo/$F- 

In a general situation the swapping time t s is deter- 
mined by the relation 2xo sm(0.5<5Ft s )/<5p = 0.57T, imply- 
ing a quickly stabilized swapping time with increasing ra- 
tio Xo/$f', see Fig. 3(a). Meanwhile the output state will 



be\a,P)±\ 



, ae lV2 ), where iff = r]it s — 0.5iT+0.58Fts 



if the inputs are two coherent states \a) and \/3). The fi- 
delity of the output state is determined by the two ratios 
\Xo/Sf\ and |r?i/xo| (given 771 = i] 2 ). As it is shown in 
Fig. 3(b), a high-quality output state will be realized 
with the proper ratios that can be achieved by adjusting 
the system parameters. 

Decoherence effect. — One limitation resulting from the 
decoherence in the parametric loop is that the transi- 
tions have some broadening and the fields detuning A, , 6 
should be much larger than the bandwidth of the corre- 
sponding transitions T g i (i — a,b and c), i.e. A, T ga g b 
and S 3> r gc . The other limitation comes from the possi- 
ble population of the other energy levels from the ground 
state \g) during the swapping process. The radiative de- 
cay in the system occurs only after the excitation of the 
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system from the ground state. As we apply a disper- 
sive interaction, the probability for the one-photon exci- 
tation scales as P a .b = g\ 2 n i,'2^i2 ( see Eqs. (|8aj) and 
(|8b[) h and the probability of the two-photon excitation 
scales asP c = \Q 1 g ly /n~[/(A 1 S) + n2g2y/n^e lAFt /(A 2 <5)| 2 
(see Eq. (p~0|) V Therefore the effective decays in the 
system should be described by the product of the de- 
cay rate of the particular state and the probability of 
the electron excitation at this level. In order to have 
the negligible contribution from the decay processes, the 
swapping time should be much shorter than the effec- 
tive radiative lifetime t s -C 1/(7^). For example, 
if one takes the swapping time for the matched four- 
photon detuning giving §p = 0, there should be the 
condition 7 Q b < - ^S 2 - for neglecting the de- 

cay of the intermediate states \a) (\b)) and the condition 
7c < f \y/nrfn+ ^n 2 e iAFt m- 1 \~ 2 S, where m = §^37 > 
for neglecting the decay of the state |c). For the example 
using 43 Ca+ with j a , b ~ 10 MHz and j c ~ 10 Hz the 



above conditions can be easily established. These condi- 
tions guarantee the outputs of the swapping process to 
be pure states. 

Summary. — With an example, we have illustrated how 
to entangle two arbitrary states l^i) and to the 

symmetric form =1= 1^2, ^l) with an induced con- 

ditional BS type coupling that avoids the physical lim- 
itation on nonlinear couplings. In contrast to all previ- 
ous schemes, the entanglement strategy is independent 
of the specific input states, e.g. SECSs and NOON states 
can be generated in the same way. The FWM process 
for realizing the effective BS coupling is within the cur- 
rent experimental technology. This approach can help to 
achieve the goals of entangling light fields with flexibility 
and creating cat states of large size. 
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